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Abstract 

We present the first measurement of two-particle momentum correlations in jets produced in 
pp collisions at -s/s = 1.96 TeV. Results are obtained for charged particles within a restricted 
cone with an opening angle of 0.5 radians around the jet axis and for events with dijet masses 
between 66 and 563 GeV/c^. A comparison of the experimental data to theoretical predictions 
obtained for partons within the framework of resummed perturbative QCD in the next-to-leading 
log approximation (NLLA) shows that the parton momentum correlations survive the hadronization 
stage of jet fragmentation, giving further support to the hypothesis of local parton-hadron duality. 
The extracted value of the NLLA parton shower cutoff scale Qeff set equal to Aqcd is found to be 
(1.4+°;?) X 100 MeV. 
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I. INTRODUCTION 



The goal of this analysis is to measure the two-particle momentum correlation in jets, 
study its dependence on jet energy, and compare the results to analytical predictions of the 
next-to-leading log approximation (NLLA) [1|. 

The evolution of jets is driven by the emission of gluons with very small transverse mo- 
menta with respect to the jet axis, i.e. less than 1 GeV/c. The theoretical predictions. 



iiA calculations 
3. NLLA pro- 



which are compared with the results of this measurement, are based on NL 
supplemented with the hypothesis of local parton-hadron duality (LPHD) 
vides an analytical description of parton shower formation, while LPHD states that the 
hadronization process takes place locally and, therefore, properties of partons and hadrons 
are closely related. Detailed studies of jet fragmentation allow one to better understand 
the relative roles of perturbative parton showering and non-perturbative hadronization in 
shaping the main jet characteristics. Past experimental studies of inclusive distributions of 
particles in jets in e~^e~ j3| and j^, 5] collisions have shown good agreement with theo- 
retical predictions, suggesting that the perturbative QCD (pQCD) stage must dominate jet 
formation, and the role of the non-perturbative stage is reduced to converting final partons 
into hadrons without significantly affecting their multiplicities and momenta. This paper 
addresses the question of whether more subtle effects, such as momentum correlation, also 
survive hadronization. The predictions for the parton momentum correlations at the level 

n 

of NLLA precision were first obtained by CP. Fong and B.R. Webber in [6[ and recently 
recalculated in the modified leading log approximation (MLLA) framework by R. Perez- 
Ramos 0]. These pQCD-driven correlations extend over a large range of parton momenta 
differences and should not be confused with phenomenological Bose-Einstein correlations 8| 
present only for parton momenta differences up to 200 MeV. 

Until now, the two-particle momentum correlations were studied only by the OPAL col- 
laboration in e+e~ collisions at a center of mass energy of ~ 91 GeV [9]. Charged particles 
in the full experimentally accessible solid angle were used in the OPAL analysis. This made 
it possible for OPAL to ignore some effects of jet reconstruction, but the correlations were 
measured over a larger cone radius than can be reliably treated theoretically. Strictly speak- 
ing, the theory describes parton shower development only within a small opening angle 6'^ 
around the jet axis, so that tan 6'c ~ Oc- The OPAL measured distributions showed a pattern 
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in qualitative agreement with theoretical predictions, but the values of the parton shower 
cutoff Qeif jll extracted from the fit of the correlation distributions were inconsistent with 
the Qeff extracted from fits of the inclusive momentum distributions jlo| . 

In this paper, we report a measurement of the two-particle momentum correlation for 
charged particles in events with dijet invariant masses in the range 66-563 GeV/c^. Mo- 
mentum correlation distributions are measured for charged particles in restricted cones with 
an opening angle of 6c = 0.5 radians around the jet axis. Events were produced at the 
Tevatron collider in pp collisions at a center of mass energy of 1.96 TeV and were recorded 
by the CDF II detector. The total integrated luminosity is 385 pb^^. 

The data sample consists of events with an expected mixture of quark and gluon jets. In 
order to compare experimental results with theoretical predictions, the analysis is carried out 
in the center-of-mass system of the two jets. The data are divided into seven bins according 
to the value of dijet invariant mass, which allows us to explore the dependence of correlation 
parameters on energy. 

The data are fit to NLLA analytical functions and the value of the parton shower cutoff 



Qeff is extracted. The corre^ 



by the pythia tune A 
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ations observed in data are compared to Monte Carlo predictions 



12| and HERWIG 6.5 [l3| event generators 



II. THEORY 



A. Next-to-leading log approximation 

NLLA allows a perturbative calculation of the parton shower by keeping all terms of 
order \n^'^ (Ej^t) and a'^\ii^"'~^{Ejet) at all orders n of perturbation theory. In these terms 
as is the strong coupling constant and Ej^t is the jet energy. Most of the particles in jets 
have /ct < 1 GeV/c [5], where fc^ is the transverse momentum with respect to the jet axis. 
Therefore, in order to successfully describe jet fragmentation, a theoretical model must be 
able to handle parton emissions at such low transverse momenta scales. 

In NLLA the requirement that parton kj- > Q cutoff , for a sufficiently high cutoff scale 
Q cutoff (typically a few GeV/c), ensures that only partons in the perturbative region are 
included in the calculation. After the resummation the value of the parameter Q cutoff can be 
lowered to the value of Aqcd- Lowering the parameter Q cutoff is equivalent to including softer 
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4: 



Q=19 GeV 
Q=50 GeV 
Q=119 GeV 



Cone 9^=0.5 




FIG. 1: NLLA inclusive parton momentum distributions for Q = EjetOc = 19, 50, and 120 GeV 
and Q^ff = 230 MeV calculated according to jl4 - 

partons in the description of the model Setting Q cutoff to its lowest value, Aqcd, maximizes 
the range of applicability of the model. The phenomenological scale replacing the two initial 
parameters Q cutoff and Aqcd is conventionally called Qeff- In theoretical calculations final 
expressions for the observables of interest are usually functions of r = ln{Q/Qcff), where 
Q = EjetOc is the so-called jet hardness. These final expressions can be expanded in powers 
of T. The Fong- Webber approach 6| keeps only terms that are fully controlled within the 
precision of the calculation, i.e. neglects all terms of order a" ln^"~^(i?jei) and higher. The 
Perez-Ramos approach [t^ partially includes higher-order terms. 

The inclusive momentum distribution function of partons in jets D(^) = ^ in NLLA is 
defined in terms of the variable ^ = In(^) where x = and p is the parton momentum. 
This distribution is predicted to have a distorted Gaussian shape [l4|: 



dN 



exp 



-I - -s6 - i(2 + 1)6^ + -s6' + —16^ 
8 2 4^ ^ 6 24 . 



(1) 



where 6 = — ^0 and ^0 is the position of the maximum of the distribution. The coefficients 
0", s, and / are the width, skewness, and kurtosis of the inclusive momentum spectrum. These 
coefficients are calculated to next-to-leading order and depend on Q^ff. Overall, the function 
has three parameters to be determined experimentally: the normalization A^, Qeff, and an 
unknown higher-order correction term 0(1) in the calculation of ^0 = 0.5r + ay^ + 0(l). 
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where a is a constant that depends on the number of colors and the number of effectively 
massless quarks. The unknown term 0(1) is expected to be independent of r. The predicted 
dependence of the inclusive momentum distribution on jet hardness is shown in Fig. [TJ 

The two-parton momentum correlation function R{C,i, ^2) is defined to be the ratio of the 
two- and one-parton momentum distribution functions: 

where D{^i,^2) = dfid(2 ' "^^^ momentum distributions are normalized as follows: 
/ D{^)dC, = (n), where (n) is the average multiplicity of partons in a jet, and 
/ D{C,i,^2)d^idC,2 = {n{n — 1)) for all pairs of partons in a jet. The average multiplicity 
of partons {n) is a function of the dijet mass Mjj and the size of the opening angle 9c- For 
6c = 0.5, {n) varies from ~ 6 to ~ 12 for Mjj in the range 80-600 GeV/c^ j^. 

The Fong- Webber approximation of Eq. for the two-parton momentum correlation 
function |[6|] can be written as follows: 

i?(Aei, A^2) = ro + ri(Aei + A^) + rslA^i - A^)', (3) 

where A^ = ^ — C,o, and the parameters ro, ri, and r2 define the strength of the correlation 
and depend on the variable r = \n{Q/Qcff). Equation ([3]) is valid only for partons with ^ 
around the peak of the inclusive parton momentum distribution, in the range A,^ ~ ±1. 
The parameters ro, ri, and r2 are calculated separately for quark and gluon jets from an 
expansion in powers of 1/ ^/T using the assumption that the number of effectively massless 
quarks Nj is 3. Keeping only terms controlled by theory, the parameters are: 

r^-175-^ r"-^ r" - -^-^ U) 

r^-133-^ r^-^ - _H (5) 

where q and g superscripts denote the correlation parameters for partons in quark jets and 
gluon jets, respectively. 

The theoretical prediction of the shape of the two-parton momentum correlation distri- 
bution function is shown in Fig. [2l Along the central diagonal A^i = — A^2! the shape of 
the two-parton momentum correlation is parabolic with a maximum at A^i = A^2- Along 
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FIG. 2: The NLLA parton momentum correlation function calculated for a gluon jet, Q = 50 GeV, 
and Qeff = 230 MeV according to 

the central diagonal A,^i = A,^2, the shape is linear and increasing toward larger values of 
A^, i.e. lower momentum partons. Therefore, the obvious features of the prediction are (1) 
the correlation should be stronger for partons with equal momenta A,^i = A,^2; and (2) the 
strength of this effect should increase for lower momentum partons. 

B. Normalization 

The correlation function from Eq. ([2]) includes two effects: (1) multiplicity fluctuations 
of partons in a jet and (2) actual momentum correlations. In this analysis, we measure 
pure momentum correlations. This can be achieved by introducing one- and two-parton 
momentum distributions normalized to unity: 

D'iO = / D'im = 1, (6) 

^'(6,6) = TT^^, / = 1. (7) 

{n[n — 1)) J 

Then, the correlation function can be defined as: 
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where F(r) = is the second binomial moment. The exphcit dependence of the 



binomial moments on the energy scale for quark and gluon jets is taken from theory |15| |: 

F,(r) = 1.75--^, F,(r) = 1.33-^. (9) 

C. Quark and gluon jets 

In theory, correlation functions are calculated for quark and gluon jets separately and 
are denoted by -Dg(0 -Dg(O) respectively. Since dijet events at the Tevatron consist of 
both quark and gluon jets, in order to compare data to theory, we rewrite the formula for 
the parton momentum distributions as follows: 

DiO = fMO + i^-h)D,iO, (10) 
= + (1 - (11) 

where fg is a fraction of gluon jets in the sample. After simple algebraic transformations, it 
can be shown that the momentum correlation Eq. ([8]) is reduced to the following: 

C(Aei, A^2) = Co + ci(Aei + A^2) + C2(Aei - A^)', (12) 
where the q coefficients {i = 0, 1, 2) are: 

r. = il!— r^j ^~ fa <? (-,o\ 

' hr'Fg + (1 - QF, ^ ^ fgV^Fg + (1 - fg)F, ^ ' ^ ^ 



where r = 1^ is the ratio of average multiplicities of partons in gluon and quark jets. The 
value of r enters in the derivation of Eqs. (jlj) and (E]) 0, Eq. (ED Q, and Eq. In 
NLLA this ratio is equal to 9/4. Values of r by the pythia 6.115 and herwig 5.6 Monte 
Carlo generators are given in 16|. 
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D. Local parton-hadron duality 



LPHD is a hadronization conjecture that suggests that the properties of hadrons and 
partons are closely related. In the simplest interpretation of LPHD, each parton at the end 
of the pQCD shower development picks up a color-matching partner from the vacuum sea 
and forms a hadron. Within LPHD, the momentum distributions of hadrons are related to 
those of partons via an energy- independent constant K^phd'- 

dNhadrons -rr partons i -\ a\ 

^^ = A-.p™-^^. (14) 

Past studies of inclusive particle distributions at e^e~ experiments {sj] and CDF 
have given strong support to the LPHD hypothesis. In this analysis, we extend the LPHD 
test by examining whether the two-particle momentum correlations predicted in the pQCD 
framework also survive the hadronization. Note that in the two-particle momentum corre- 
lation given by Eq. ([2]), K^phd simply cancels, suggesting that the correlation distributions 
for hadrons and partons are expected to be the same. 

III. CDF II DETECTOR 

Data used in this analysis were recorded with the CDF II detector. The detector was 
designed for precision measurements of the energy, momentum and position of particles 
produced in proton-antiproton collisions. This section provides a brief overview of the com- 



ponents re 
elsewhere 



evant to our analysis. A detailed description of the entire detector can be found 



17|. 



CDF II uses a cylindrical coordinate system with the positive z direction selected along 
the proton beam direction and azimuthal angle measured around the beam axis. The 
polar angle 6 is measured with respect to the positive z direction and the pseudorapidity r] 
is defined as ?7 = — In tan(|) . 

The CDF II tracking system is placed inside a 1.4 T solenoidal magnet. A Layer 00 
single-sided silicon microstrip detector is mounted directly on the beam pipe, at an inner 
radius of 1.15 cm and an outer radius of 2.1 cm. A five-layer silicon microstrip detector 
(SVX II) is closest to the beamline, and is situated at a radial distance of 2.5 to 11 cm 
from the beam. The SVX II consists of three separate barrel modules with a combined 
length of 96 cm. Three of the five layers combine a r-(j) measurement with a 2;-coordinate 
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measurement while the remaining two layers combine r-0 with a small stereo angle of 1.2°. 
Three additional intermediate silicon layers (ISL) are positioned between 19 and 30 cm. 
The SVX II is surrounded by the central outer tracker (COT), an open-cell drift chamber 
providing up to 96 measurements of a charged particle track over the radial region from 40 
to 137 cm. The 96 COT measurements come from 8 superlayers of 12 sense wires each. The 
superlayers alternate between axial and 3° stereo. The pseudorapidity region covered by the 
COT is |?7| < 1.0. The momentum of a charged particle is determined by the curvature of 
its trajectory in the magnetic field. 

The CDF II tracking system is surrounded by calorimeters used to measure the energy 
of charged and neutral particles. In the central region the central electromagnetic (CEM), 
central hadronic (CHA), and wall hadronic calorimeters are made of lead (electromagnetic) 
and iron (hadronic) layers interspersed with scintillator. The pseudorapidity region covered 
by these calorimeters is |?7| < 1.3. The segmentation of the central calorimeters is 15° 
in (j) and 0.1 units in rj. The measured energy resolutions for the CEM and CHA are 
(y{E)/E = 13.5%//E^©2% and a{E)/E = 75%/v^©3%, respectively. Here Et = EsinO 
is transverse energy of an incident particle (electron for CEM and pion for CHA) and energies 
are measured in GeV. 



IV. EVENT SELECTION 



A. Triggers 

Events were collected using a single-tower trigger [l^ with a transverse energy Et thresh- 
old of 5 GeV (ST05) and with single jet triggers with Et thresholds of 20 (J020), 50 (J050), 
70 (J070), and 100 (JlOO) GeV. Each of the jet triggers had a different sampling rate so as 
not to saturate the available trigger bandwidth. 



B. Jet reconstruction algorithm 



Jets are reconstructed based on the calorimeter information using a cone algorithm |19 |. 
The algorithm starts with the highest transverse energy tower and forms preclusters from 
an unbroken chain of continuous seed towers with transverse energy above 1 GeV within a 
window of 7 X 7 towers centered at the originating seed tower. If a seed tower is outside 
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this window, it is used to form a new precluster. The coordinates of each precluster are the 
ii^T-weighted sums of (f) and t] of the seed towers within this precluster. In the next step, all 



towers with Et > 0.1 GeV within R = ^ (A0)^ + (Ar/)^ = 1.0 of the precluster are merged 
into a cluster, and its [r], 0)-coordinates are recalculated. This procedure of calculating 
cluster coordinates is iterated until a stable set of clusters is obtained. A cluster is stable 
when the tower list is unchanged from one iteration to the next. If the clusters have some 
finite overlap, then an overlap fraction is computed as the sum of the Et of the common 
towers divided by the Et of the smaller cluster. If the fraction is above a cutoff value equal 
to 0.75, then the two clusters are combined. If the fraction is less than the cutoff, the shared 
towers are assigned to the closest cluster. The raw energy of a jet is the sum of the energies 
of the towers belonging to the corresponding cluster. Corrections are applied to the raw 
energy to compensate for the non-linearity and non-uniformity of the energy response of the 
calorimeter, the energy deposited inside the jet cone from sources other than the leading 
parton, and the leading parton energy deposited outside the jet cone. A detailed description 



of this procedure can be found in |20|]. 



C. Offline selection 

Cosmic ray events are rejected by applying a cutoff on the missing transverse energy J^t 



significance 2l(|, defined as ]^t / V^Et, where T,Et = ^{Et is the total transverse energy 
of the event, as measured using calorimeter towers with E^ above 100 MeV. The thresholds 
are 3.0 GeV^^^ for data collected using a single tower trigger with Et threshold of 5 GeV, 
and 3.5, 5.0, 6.0, and 7.0 GeV^/^ for data collected using jet triggers with thresholds of 20, 
50, 70, and 100 GeV, respectively. 

To ensure fully efficient vertex and track reconstruction, we require only one vertex in 
the event with \z\ < 60 cm. 

To ensure robust and high efficiency track reconstruction and applicability of the back- 
ground removal technique (see Sec. V(B)), only events with both leading jets in the central 
region (Ir^l < 0.9) are selected. To reject events with poorly measured jets, we require the 
two leading jets to be well balanced in Et'- k±/ {E^ + E^) < 2ak^. Here E\, and E'^ are the 
transverse energies of the first and the second leading jets, respectively, and fc_L is: 



ki. = - kl (15) 
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FIG. 3: Definition of variables for the jet balance requirement. The plane perpendicular to the 
beamline is shown. The vector k represents a vector sum of the two leading jets' momenta. The k^^ 
and k± components of k are parallel and perpendicular to the bisector of the two jets, respectively. 

kii = {E^ + E^) ■cos{tlj/2), (16) 

where k is a vector sum of momenta of the two leading jets, ip is the angle between the two 
leading jets, and akj^ is the resolution of kj_. The definitions of k, and k\\ are illustrated 
in Fig. [3l The component k± is known to be sensitive to jet energy mismeasurements, while 
is more sensitive to hard gluon radiation. 

In events with high energy jets, a single particle emerging from a jet at a sufficiently large 
angle with respect to the jet axis can be mistakenly identified as a separate jet. A jet can 
also be produced from the underlying event. Therefore, rejection of all events with more 
than two jets can introduce possible biases. We allow up to two extra jets, but their energy 
is required to be small: < 5.5 GeV +0.065(-E^ + E^), where is the transverse 

energy of an extra jet. 

After application of the event selection cuts, the final sample consists of approximately 
250,000 events and is further divided into seven bins according to the dijet mass as measured 
by the calorimeters and defined as: 

M,, = ^/{E^+E2}yJ^{K+W/c', (17) 

where E and P are the energies and momenta of the two leading jets, respectively. 

The mass bin boundaries, average invariant mass (Mjj), and number of events in each bin 
are given in Table [B The bin width is selected to be 3 ■ ctm,, , where aM,, is the calorimeter 
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TABLE I: Dijet mass bins boundaries, average dijet invariant mass {Mjj), average E'jet-weighted 
jet hardness Q, and number of events in each bin after the event selection requirements Ng^. 



Bin Trigger Low edge (GeV/c^) High edge (GeV/c^) {Mjj) (GeV/c^) Q (GeV) N 



ev 



1 ST05 66 95 76 19 15229 

2 J020 95 132 108 27 77246 

3 J020 132 180 149 37 17682 

4 J050 180 243 202 50 80608 

5 J050 243 323 272 68 18528 

6 J070 323 428 361 90 12000 

7 JlOO 428 563 475 119 19150 



resolution for the dijet mass determination, ~10-15%. Events with dijet invariant mass 
below 66 GeV/c^ are not used in the measurement because of the lower trigger efficiency. 



D. Systematic uncertainties associated with the event selection 

The sensitivity of the two-particle momentum correlation parameters Cq, ci, and C2 to var- 
ious uncertainties in the event selection procedure is evaluated as follows. For each source of 
systematic uncertainty the so-called "default" and "deviated" two-particle momentum cor- 
relation distributions are obtained. The default distribution is produced using the standard 
selection requirements described in this article. Then, the deviated distribution is obtained 
by varying all relevant parameters according to the estimated systematic uncertainty (one 
source of uncertainty at a time). For each bin in correlation C(A,^i, A{2), a scale factor is 
calculated by taking the bin-by-bin ratio of the deviated and default distributions: 

^ _ C{A^i, ^£,2) deviated ^^^-^ 
C(A^1, ^^2) default 

The difference between correlation distributions in the data with and without this bin- 
by-bin scale factor is taken as a measure of the systematic uncertainty: 



AC(Aa, Ai2)Data = |(1 - e) " C(Aa, Ai2)Data 
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(19) 



TABLE II: Summary of the systematic uncertainties of the correlation parameters cq, ci, and C2 
for the dijet mass bin with Q = 50 GeV. 



Origin of systematic uncertainty 


Acq 


Aci 


AC2 


Luminosity dependence 


0.001 


0.004 


0.002 


Jet energy scale 


0.001 


0.001 


0.001 


Balance and extra jet cuts 


0.006 


0.001 


0.003 


Mismeas. of jet direction 


0.006 


0.008 


0.007 


Track quality cuts 


0.014 


0.008 


0.006 


Underlying event background 


0.001 


0.004 


0.001 


Tracking inefficiency 


0.011 


0.001 


0.002 


Neutral particles 


0.002 


0.002 


0.001 



Further in this section, we discuss sources of systematic uncertainties at the level of the 
event selection; their contributions to the values of cq, ci, and C2 are given in Table HTl 

In each trigger sample only the events with trigger efficiency higher than 99% are used. 
To check that trigger effects do not bias the measurement, we verify the continuity of the 
distributions of particle multiplicity in a jet in the transition between adjacent dijet trigger 
samples. No detectable offsets are observed. 

To evaluate the uncertainty due to the value of the parameter R of the jet reconstruction 
algorithm, we compare the results of the measurement using three different values of R 
(0.4, 0.7, 1.0). This effect proved to be small compared to the other sources of systematic 
uncertainty. 

We require only one vertex in the event, which greatly reduces the contribution of multiple 
proton-antiproton interactions in the same bunch crossing. However, in some cases two 
vertices can be very close to each other and be reconstructed as a single vertex. This 
can become significant at high instantaneous luminosity. To evaluate the uncertainty due 
to this effect, we divide each dijet mass bin into smaller bins based on the instantaneous 
luminosity. Momentum correlation distributions are compared in these smaller samples and 
the difference is taken as a measure of the systematic uncertainty. 

To evaluate the uncertainty due to the jet energy corrections we use parameterizations in 
which the jet energy scale is shifted by plus or minus one standard deviation [2^. We then 
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reclassify the events according to their dijet mass. The difference between the default and 
deviated distribution is assigned to be the systematic uncertainty. 

We use Monte Carlo dijet samples produced by pythia tune A to study systematic 
uncertainties associated with the jet balance requirement, the number of allowed extra jets, 
and their energy. The default two-particle momentum correlation distribution is compared 
to the deviated one. The deviated distribution has no requirements imposed on the jet 
balance or on the extra jet number and extra jet energy. 

Monte Carlo simulations are used to evaluate the systematic uncertainty due to mismea- 
surement of the jet direction. Two-particle momentum correlations are compared for two 
cases. In one case particles are counted in a restricted cone around the jet direction as 
determined by the detector response in the simulation. In the second case the direction of 
primary partons from the hard scattering as given by pythia tune A is used for the cone 
axis. 

V. TRACK SELECTION, CORRECTIONS, SYSTEMATIC UNCERTAINTIES 

Measurements described below are performed in the dijet center of mass frame. For 
Lorentz boosts all particles arc treated as pions. Experimentally wc define the variable ^ 
as £ = Infl/x) = In-^^, where Ej^t is the jet energy as measured by the calorimeters 

Ptrack 

and Ptrack is the track momentum as measured by the tracking system. The correlation 
distributions are measured for all track pairs that pass track quality requirements and lie 
within a restricted cone of opening angle 9c — 0.5 radians relative to the jet axis. The peak 
position of the inclusive momentum distribution is constant for a given jet hardness Q 

and is obtained from the data. The measurements are corrected for various backgrounds 
both correlated and uncorrelated with jet direction. 

A. Track quality requirements 

Several selection requirements are applied to ensure that the tracks in the measurement 
originate at the primary vertex and are not produced by cosmic rays, multiple pp interac- 
tions within the same bunch crossing, 7-conversions, and A-decays, or other types of 
backgrounds. 
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In our analysis we require full three-dimensional track reconstruction. The description 



of CDF II track reconstruction can be found in 
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23| . Poorly reconstructed and spurious 
tracks are removed by requiring a good track fitting parameter Xcot < 6.0. Charged 
particles are required to have transverse momentum px > 0.3 GeV/c. 

The parameter Az is defined as the difference between the z position of the track at the 
point of its closest approach to the beamline and the z position of the primary vertex. This 
parameter is used to remove tracks not originating at the primary interaction by requiring 
\Az\ < S-cta^j where a^z is determined for different categories of tracks based on the number 
of SVX II, ISL, and COT hits. 

Tracks produced from 7-conversions are removed using a combination of requirements 
on impact parameter do and the distance Rconv (see Fig. H]). The impact parameter do is 
defined as the shortest distance in the r — (p plane between the beamline and the trajectory 
of the particle obtained by the tracking algorithm fit. It can be shown that for electrons 
and positrons originating from 7-conversion: 



«... « (20) 

where pt is the transverse momentum of the charged particle in GeV/c, B is the magnetic 
field in Tesla and Rconv is measured in meters. Monte Carlo studies indicate that the 
combined requirement of |(io| < 5 ■ a do or Rconv < 0.13 m is more efficient at removing 
7-conversion tracks than the do requirement alone. The resolution of the impact parameter, 
ado, is evaluated for different categories of tracks based on the number of SVX II, ISL, 
and COT hits. The value Rconv = 0.13 m is motivated by the location of SVX II readout 
electronics. Indeed, conversions occurring at this radius are clearly seen in the data. 

To verify the effectiveness of the track quality cuts, we compare distributions of the 
inclusive particle multiplicity and momentum in pythia tune A at the generator level and 
at the level of the detector simulation (cdfsim) . The comparison is shown in Figs. [5] and 
ini CDFSIM propagates particles through the detector including both conversions and in- 
flight decays to simulate the CDF II detector response. The agreement after selection cuts 
are applied confirms that the cuts do remove most of the background tracks. The effect 
of the remaining fraction of secondary tracks is estimated by comparing the correlation 
distributions C(A^i, A^2) at the charged hadron level and the cdfsim level and producing 
a corresponding bin-by-bin scale factor. The difference between distributions in data with 
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beam line 

FIG. 4: Schematic illustration of the distance Rconv from the beamline to the point where the 
conversion occurred. Here, do is the impact parameter. 
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FIG. 5: Monte Carlo track multiplicity in jets before and after applying track quality cuts. The 
distributions are for the dijet mass bin with Q = 50 GeV. Particles are counted within a cone of 
opening angle 6c = 0.5 radians. CDFSIM refers to the full CDF II detector simulation. 

and without this scale factor applied is assigned as the systematic uncertainty associated 
with the track quality cuts. 



B. Underlying event background subtraction 

Generally, tracks from the underlying event tend to dilute the two-particle momentum 
correlation. It is not possible to correct for this effect on an event-by-event basis, but 
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FIG. 6: Inclusive momentum distributions of Monte Carlo tracks in jets before and after applying 
track quality cuts. The distributions are for the dijet mass bin with Q = 50 GeV. Particles are 
counted within a cone of opening angle 9c = 0.5 radians, cdfsim refers to the full CDF II detector 
simulation. 

the average correction factor can be reconstructed statistically. In order to correct for the 
underlying event contribution, we apply the following procedure. On an event-by-event 
basis, two complementary cones are positioned at the same polar angle with respect to the 
beamline as the original dijet axis but in the plane perpendicular to the dijet axis as shown 
in Fig. [71 Complementary cones defined this way are at 90° in (i.e. as far as possible) 
from the dijet axis. This can be done when the dijet axis is within 45° < 9 < 135°, and this 
condition is automatically satisfied by our event selection. We assume that cones formed in 
such a fashion collect statistically the same amount of background (which is uncorrelated 
with jets) as the cones around the jet axis j^l. 

In order to obtain the corrected expression for C(A,^i, A,^2), one needs to subtract the 
background from the one- and two-particle momentum distributions. This can be achieved 
by considering particles in jet cones together with particles in complementary cones. It can 
be shown that the momentum distributions after background subtraction D are: 



D(0=/^,et(0-^co„.pKO, (21) 
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FIG. 7: Illustration of the definition of complementary cones. The unlabeled arrows are the axes 
of the cones complementary to jets 1 and 2. The complementary cone makes the same angle 9 with 
the beamline as the jet cone. 

where the jet subscript denotes the distribution for particles in jet cones, compl denotes the 
distribution for particles in complementary cones, and jet + compl denotes the distribution 
for the combined set of particles in either jet cones or complementary cones. 

To evaluate systematic uncertainties associated with the background subtraction using 
the complementary cones, we use the following procedure. The amount of background in 
a jet cone is increased by a factor of two by adding tracks from the complementary cone 
of another event. Then, the background subtraction procedure described above is applied 
taking into account the artificially doubled background. After the subtraction the correlation 
distribution is expected to be the same as the distribution using the original background. 
The difference between the two-particle momentum correlation distributions obtained after 
the subtraction of either the original or the doubled background is assigned as a measure of 
the systematic uncertainty. 

C. Tracking inefficiency 

A high efficiency of track reconstruction is ensured by selecting events with central jets. 
However, there still may be non-reconstructed tracks inside the jet. To evaluate the corre- 
sponding systematic uncertainty, we have modeled the track reconstruction inefficiency using 
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the function P{C) = pi + P2^-, which denotes the probabihty of losing a track with given ^. 
Values of the parameters pi and p2 were varied over a range far exceeding the estimated 
COT inefficiency. The correlation distributions show a very weak dependence on tracking 
inefficiency. The range of momentum correlation variation in this tracking inefficiency model 
is taken as a measure of the systematic uncertainty (see Table HTj) . 

D. Neutral particles 

Theoretical predictions of correlation distributions are done at the parton level, while 
LPHD relates final partons to hadrons, assuming that all hadrons are counted. The analysis, 
however, is done for charged particles only. To estimate the effect of neutral particles the 
momentum correlation in a PYTHIA tune A sample is compared for charged particles and all 
particles. The difference is assigned as the corresponding systematic uncertainty (see Table 

m- 

E. Resonance decays 

The presence of resonance decays may be expected to cause differences between the 
correlation in data and the theoretical predictions. We examine this effect by comparing 
the correlations in Monte Carlo events for hadrons before and after resonance decays. We 
find that this results in insignificant changes in C(A^i, A^2) and does not change the overall 
level of the correlation. 

F. Heavy flavor jets 

Theoretical predictions of correlation distributions are obtained for jets originated from 
gluons or light quarks only. In the data sample we expect a small fraction (~ 5%) of heavy 
fiavor jets. To estimate the size of this effect we repeat the analysis with the assumption that 
the correlations in heavy fiavor jets are same as in gluon jets. This translates into a 3 MeV 
change in the value of Qejj and is negligibly small compared to the size of the systematic 
uncertainty. 
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VI. NLLA FITS TO DATA 



The inclusive momentum distributions D(^) = ^ in all seven experimental dijet mass 
bins are simultaneously fit to the theoretical Fong-Webber function. In the fit the Qeff 
and 0(1) parameters are required to have same value in all dijet mass bins while the nor- 
malization parameter N{t) is allowed to vary from one bin to another. Figure [8] shows the 
distributions in data corresponding to the dijet mass bins with Q = 27, 50, and 90 GeV. The 
error bars correspond to both the statistical and systematic uncertainties added in quadra- 
ture. The solid curves correspond to the fit of the data to the theoretical Fong-Webber 
function, and the dashed curves represent the extrapolations out of the fit regions. The 
extracted values of the fit parameters are Qejj = 180 ± 40 MeV and 0(1) = —0.6 ± 0.1, 
where the uncertainties are statistical and systematic added in quadrature. The value of 
Qeff is consistent with the results of a previous CDF measurement Sj. 

The two-particle momentum correlation distributions 0(A^i, A^2) are produced for seven 
bins of dijet mass and do show the shape predicted by theory. In this paper we plot the 
central diagonal profiles A^i = — A^2 and A^i = A^2 (see Fig. [21) of the distributions. 
Figures [HI [lOl and [H] show the distributions corresponding to the dijet mass bins with 
Q = 27, 50, and 90 GeV, respectively. The bin size A^ = 0.2 is chosen to be much wider 
than the momentum resolution in the fitted range. The smaller error bars correspond to 
the statistical uncertainty only, while the larger error bars correspond to both the statistical 
and systematic uncertainties added in quadrature. The 2-dimensional momentum correlation 
distribution is fit according to Eq. (|T2l) with three free parameters cq, ci, and C2- The solid 
lines in Figs. M, [lOl and [H] show the profiles of the fit functions. The extracted values of the 
fit parameters are given in Table [Till The fit range — 1 < A^ < 1 is motivated by the region 
of validity of the NLLA calculations. 

The dash-dotted lines in Figs. [9], [TOl and [11] correspond to the theoretical curves given 
by Eq. ( JT2l) for Q^ff = 180 ± 40 MeV, extracted from fits of the inclusive momentum distri- 
butions. The dashed lines correspond to the results of the Perez-Ramos calculation for the 
value of Oeff = 230 ± 40 MeV extracted from fits of the inclusive momentum distributions 

n 

to the MLLA function [5|. The fraction of gluon jets in the sample, used to model the theo- 
retical prediction for quark and gluon jets, is obtained using pythia tune A with CTEQ5L 
parton distribution functions 26|. 
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TABLE III: Summary of the correlation parameters cq, ci, and C2 measured in seven dijet mass 
bins. The first uncertainty is statistical, the second one is systematic. 



Q (GeV) 


Co 


Cl 


C2 


19 


1.078 ±0.007 ±0.016 


0.081 ± 0.006 ± 0.016 


-0.047 ± 0.006 ± 0.008 


27 


1.076 ± 0.003 ± 0.022 


0.068 ± 0.002 ± 0.015 


-0.038 ± 0.002 ± 0.012 


37 


1.075 ±0.005 ±0.018 


0.057 ±0.004 ±0.013 


-0.031 ±0.004 ±0.012 


50 


1.079 ±0.002 ±0.019 


0.051 ±0.002 ±0.014 


-0.029 ± 0.002 ± 0.010 


68 


1.081 ± 0.004 ± 0.028 


0.040 ± 0.004 ± 0.012 


-0.027 ±0.004 ±0.011 


90 


1.081 ± 0.005 ± 0.023 


0.046 ± 0.004 ± 0.015 


-0.024 ± 0.004 ± 0.014 


119 


1.077 ± 0.004 ± 0.033 


0.028 ±0.003 ±0.013 


-0.019 ±0.003 ±0.015 



The systematic uncertainty due to the parton distribution functions is evaluated by com^ 



paring results for the fraction of gluon jets fg obtained using CTEQ5L and CTEQ6.1 [27 1 
PDF sets. The systematic uncertainty due to the value of r was evaluated by taking the 
difference between the theoretical value {r theory = 9/4), used as default, and experimental 
value {vexp = 1-8) [3] and propagating it to the value of Qeff- Both systematic uncertainties 
were found to be negligible. 

The overall qualitative agreement between the data and the Fong- Webber calculation 61] 
is very good. The data follow the theoretical trends and show an enhanced probability of 
finding two particles with the same value of momenta (indicated by the parabolic shape 
of the A,^i = — A^2 central diagonal profile with its maximum at A,^i = A,^2 = 0). This 
effect becomes larger for particles with lower momenta (the positive slope of the A^i = A^2 
central diagonal profile). An offset in the overall level of correlation is observed in all seven 
dijet mass bins, indicating that the Fong- Webber prediction overestimates the parameter cq 
of the correlation. The Perez-Ramos curves [7] qualitatively show the same trends; however, 
the quantitative disagreement is obviously larger for the Perez- Ramos predictions compared 
to the Fong- Webber predictions 

Figure [T2l shows the dependence of parameters cq, ci, and C2 on jet hardness Q. Each data 
point corresponds to the value of one parameter measured in a particular dijet mass bin. The 
Co parameter shows almost no dependence on Q, while |ci| and |c2| decrease with increasing 
Q. This indicates that the correlations are stronger in low energy jets. The distributions 
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FIG. 8: Inclusive momentum distributions of particles in jets in the restricted cone of size 9c = 0.5 
radians for the dijet mass bin with Q = 21 GeV (top), Q = 50 GeV (middle), and Q = 90 GeV 
(bottom). The solid curves correspond to the fit of CDF data to the theoretical Fong- Webber 
function (as calculated in |^]), the dashed curves represent the extrapolations out of the fit regions. 
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FIG. 9: Two-particle momentum correlations in jets in the restricted cone of size 9c = 0.5 radians 
for the dijet mass bin with Q = 21 GeV (top). Central diagonal profiles A^i = — A,^2 (middle) and 
A^i = A,^2 (bottom) of the distributions are shown. The correlation in data is compared to that 
of theory (as calculated in [6] for Q^g = 180 MeV and in P] for Q^s = 230 MeV). 
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FIG. 10: Same as in Fig. M but for Q = 50 GeV. 
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FIG. 11: Same as in Fig. [9] but for Q = 90 GeV. 
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FIG. 12: The dependence of correlation parameters C2 (top), ci (middle), and cq (bottom) on 

jet hardness. The Fong- Webber function is fit to the CDF data points. The Fong-Webber 

predictions for pure quark and pure gluon jet samples are also shown. 
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are fit to the Fong- Webber function witfi Qeff treated as tlie only free parameter. The 
fits are represented by solid lines. Theoretical curves for pure quark and gluon jets in the 
final state are also shown. We use the results of the Fong- Webber calculation to fit 
the dependence of these parameters on jet hardness and to extract the parameter Qejj- 
Results of the Perez- Ramos calculation are not used for the measurement of Qejj due to 
the lack of the corresponding analytical expressions. The value of Qejf obtained from the 
fit of Ci is 145 ± 10(stat)^g5(syst) MeV. The value of Qejj obtained from the fit of C2 is 
129 ± 12(stat)+7J(syst) MeV. The average value of Qeff extracted from the combined fit 
of Ci and C2 is IST^gg MeV and is consistent with Qeff extracted from the fits of inclusive 
particle momentum distributions. The dependence of Cq on Q has an offset of ~0.06. This 
parameter, as opposed to Ci and C2, is very sensitive to the peak position of the inclusive 
momentum distribution. In the data the correlation distributions are measured around the 
true peak position while in the theoretical calculation of the unknown constant term 0(1) 
as well as all terms beyond the leading order are neglected. Therefore, theory can control 
only the dependence of this parameter on energy and not its absolute value. For this reason 
we exclude Cq from the measurement of Qeff- A formal fit of the dependence of Cq on Q to 
the theoretical function gives the value Qeff = 0.10 ± 0.08 MeV. This value, however, does 
not have physical meaning for the above mentioned reasons. Other than the offset, cq shows 
very weak, if any, Q dependence, which is consistent with the theory. As a cross-check we 
have measured correlation distributions for pairs of tracks from opposite jets. For our value 
of the opening angle 6c = 0.5 radians, no correlations are observed. 



VII. COMPARISON TO MONTE CARLO 



We compare the momentum correlation distributions of charged particles in data to 
PYTHIA tune A and HERWIG 6.5 predictions. Predictions of the two Monte Carlo generators 
are in good agreement with each other and with results obtained from data. Figures [131 
[TH and [15] show the correlation distributions in data compared to pythia tune A and 
HERWIG 6.5 predictions at the level of stable charged hadrons. 
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FIG. 13: Hadron-level two-particle momentum correlations in jets in the restricted cone of size 
9c = 0.5 radians for the dijet mass bin with Q = 21 GeV using pythia tune A (top). Data 
correlations are compared to the hadron momentum correlations using the PYTHIA tune A and 
HERWIG 6.5 event generators. Central diagonal profiles A^i = — A^2 (middle) and A^i = A^2 
(bottom) of the distributions are shown. 34 
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FIG. 14: Same as in Fig. [T3 but for Q = 50 GeV. 
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FIG. 15: Same as in Fig. US but for Q = 90 GeV. 
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VIII. SUMMARY 



The two-particle momentum correlation distributions of charged particles in jets from 
dijet events have been measured over a wide range of dijet masses from 66 to 563 GeV/c^. 
The jets were produced in pp collisions at a center of mass energy of 1.96 TeV. The mea- 
surements have been performed for particles in a restricted cone around the jet direction 
with opening angle 6c = 0.5 radians. 

The data are compared to the next-to-leading log approximation calculations combined 
with the hypothesis of local parton-hadron duality (LPHD). Overall, the data and the theory 
show the same trends over the entire range of dijet energies. The parton shower cutoff scale 
Qeff is set equal to Aqqd and is extracted from fits of the dependence of the correlation 
parameters, ci and C2, defining the strength of the correlation, on jet hardness Q. The 
average value of Qejj extracted from the combined fit of Ci and C2 is 137tg9 MeV and is 
consistent with Qgff extracted from the fits of inclusive particle momentum distributions 
and with the results of a previous CDF measurement js;]. As predicted, the parameter cq 
shows little, if any, dependence on jet energy; however, we observe a substantial systematic 
offset between the experimental and theoretical values. The parameter cq is excluded from 
the measurement of Qeff because of its large theoretical uncertainty. The modified leading 
log approximation predictions qualitatively show the same trends; however, the quantitative 
disagreement with the data is obviously larger in this case. 

The PYTHIA tune A and herwig 6.5 Monte Carlo event generators are found to reproduce 
the correlations in data fairly well. 

The results of this analysis indicate that the parton momentum correlations do survive 
the hadronization stage of jet fragmentation, giving further support to the hypothesis of 
LPHD. 
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